Phase diagram for ultracold bosons in optical lattices and superlattices 
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We present an analytic description of the finite-temperature phase diagram of the Bose-Hubbard 
model, successfully describing the physics of cold bosonic atoms trapped in optical lattices and su- 
perlattices. Based on a standard statistical mechanics approach, we provide the exact expression for 
the boundary between the superfluid and the normal fluid by solving the self-consistency equations 
involved in the mean-field approximation to the Bose-Hubbard model. The zero-temperature limit 
of such result supplies an analytic expression for the Mott lobes of superlattices, characterized by a 
critical fractional filling. 
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I. INTRODUCTION 

The standing wave produced by the interference among 
counterpropagating laser beams gives rise to a periodic 
potential commonly used to fragment and trap clouds 
of ultracold (possibly condensed) alkali atoms 0, 0, Q • 
The local minima of such trapping potential are the sites 
of the so-called optical lattice p|. The use of multiple 
wavelength laser beams allows to obtain superlattices, 
namely more structured periodic potentials characterized 
by a spatial modulation of the depth of the lattice wells 

In the case of bosonic atoms cooled to within the lowest 
Bloch band of the periodic potential, it can be shown Q| 
that the physics of the system is described by the well- 
known Bose-Hubbard Hamiltonian |lfj| 
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where the subscripts j, h label the sites of the optical 
lattice, (aj) creates (annihilates) a boson at site j, 
rij = djCLj counts the bosons at site and the symbol 
~ restricts the sum over h to the nearest neighbours of 
j. This is obtained by a scheme analogous to the tight- 
binding approximation commonly adopted for the study 
of electrons in solids, i.e. by expanding the state of the 
system onto a set of wavefunctions localized at the local 
minima of the trapping potential. The parameters U, t, 
Tjh and Vj are hence given in terms of overlap integrals 
between the localized wavefunctions at neighbouring sites 
and the trapping potential. Specifically, U represents the 
repulsive boson-boson interaction, Vj is the local poten- 
tial at site j and t ■ Tjh is the hopping amplitude between 
sites j and h. More in detail, t is a global scaling fac- 
tor determined by the laser intensity, while Tjh = Thj is 
a local parameter depending on the details of the opti- 
cal potential in the region between lattice sites j and h 
(see Fig. [3] for a schematic representation of the optical 



superlattices considered in the following). A fine tuning 
of parameters Vj, t and Tij can be in principle obtained 
via suitable variations of experimental parameters such 
as the intensity, the frequency and the geometric setup 
of the laser beams producing the optical (super)lattice 
For generic superlattices Tjh and Vj are periodic 
functions of the lattice labels , whereas in the case of 
the usual single-wavelength optical lattices they are usu- 
ally assumed to be site-independent, so that one can set 
Tjh = 1 and Vj = without loss of generality. 

The parameter /i appearing in Hamiltonian Q is the 
usual chemical potential of the grand canonical statistical 
approach, and is fixed by the total number of bosons in 
the system. 

Hamiltonian is also strictly related to systems other 
than the one under concern, such as Josephson junction 
arrays and quantum spin systems on lattices |12lll3l |. The 
hallmark of such class of systems is no doubt the quan- 
tum phase transition between a superfluid and a (Mott) 
insulator phase E3 originating from the competition be- 
tween the repulsive and kinetic term of the Hamiltonian, 
whose magnitude are proportional to the parameters U 
and t, respectively. The fine tuning of these parameters 
made possible by the striking progress in optical lattice 
techniques allowed Greiner and co-workers to observe the 
superfluid-insulator transition in a recent breakthrough 
experiment |14| . More in general, ultracold neutral atoms 
make an ideal benchmark for testing the properties of 
widespread models of condensed matter physics EH . 

It is worth recalling that the above quantum phase 
transition is rigorously present only at zero temperature 
[16|, whereas at finite temperature thermal fluctuations 
induce a classical phase transition between a superfluid 
and a normal phase. However, at sufficiently low tem- 
peratures, a remnant of the insulating phase still persists 
within the normal phase. Indeed in these conditions it is 
possible to observe a sharp crossover between a compress- 
ible normal fluid and a phase characterized by a vanishing 
compressibility, which, for all practical purposes, can be 
considered a Mott insulator 
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The zero-temperature phase diagram of the Bose- 
Hubbard model has been widely studied using a vari- 
ety of techniques, including the mean-field approaches 
EE El IIS |21| . strong coupling perturbative expansion 
[221 l'2.'l |2 1| , density matrix renormalization group l25l 
and of course Quantum Monte Carlo simulations [2fiLl27l |. 

As to the finite-temperature case, some early numerical 
results concerning the homogeneous lattice are reported 
in Ref. |2^, where a coarse graining mean-field approach 
is adopted, and in Ref. ITtI based on a random phase ap- 
proximation refining the mean-field approach therein pro- 
posed. In Ref. [22l the weak-repulsion limit is addressed, 
and some results are obtained numerically within a dif- 
ferent mean-field scheme, based on the linearization of 
the repulsive term rather than on the decoupling of the 
hopping term, like in Ref. El Quite recently, Dicker- 
scheid and co-workers 0] adopted a slave-boson tech- 
nique allowing to include the finite temperature effects 
in the mean-field picture of Ref. ljj The ongoing inter- 
est in the issue under examination is further confirmed 
by Ref. 13(3 — where a multiband model is addressed — 
and Ref. |3l| — where some analytic results are obtained 
by interpolating two different perturbative schemes and 
subsequently checked against density matrix renormal- 
ization group simulations. All of the above listed results 
have been obtained either numerically or applying some 
further approximation such as introducing tight restric- 
tions on the number of particles per site. Note however 
that, despite such restrictions, the latter approach may 
prove sufficient to give satisfactory results within circum- 
scribed regions of the phase diagram. 

Here we focus on the mean-field approach to Hamilto- 
nian 0) proposed by Sheshadri et al. 0|j anci ) f° r an y 
temperature T, we determine analytically the boundary 
of the superfluid domain of the phase diagram thereof. 
In this framework the phases of the system are char- 
acterized in terms of the so-called superfluid order pa- 
rameter, to be determined as the stable fixed point of 
a self-consistency equation. Detailedly, this parameter 
vanishes in the normal fluid phase, whereas it has a fi- 
nite value in the superfluid phase. We determine the 
critical boundary between these phases by discussing the 
(parameter dependent) stability of the fixed point corre- 
sponding to the normal phase. Furthermore we discuss 
the above-mentioned crossover between the compressible 
normal fluid and the insulator-like phase taking place 
outside the superfluid domain. Other than the usual d- 
dimensional homogeneous lattice, we consider a generic 
one-dimensional ^-periodic superlattice, providing a so- 
lution in terms of the maximal eigenvalue of a I x l- 
matrix. Explicit results are given for the the 2-periodic 
and for a special case of the 3-periodic superlattice, where 
such maximal eigenvalue can be easily worked out. The 
zero-temperature phase-diagram of the above mentioned 
systems is recovered taking the appropriate limit in our 
results. In particular, for superlattices, we find that ra- 
tional filling lobes appear besides the usual integer-filling 
Mott domains. Also, we observe that the occurrence of 



the latter can be prevented with suitable choices of the 
supercell potential profile. All of our results prove equiv- 
alent to those obtained adopting the standard numerical 
algorithm, based on a self-consistent iterative procedure. 

The plan of this paper is as follows. In section ITU we 
briefly recall the mean-field approach presented in Ref.llTl 
and introduce the finite-temperature self-consistency 
condition. In Section IIIII we shortly address the homo- 
geneous lattice case, providing the exact expression for 
the boundary of the superfluid domain, and discussing 
the crossover between the compressible normal fluid and 
a insulator-like phase. These results are extended to the 
case of periodic superlattices in Sec. IIVI where the two 
above-mentioned special cases are explicitly considered. 
Most of the technical details of our derivation are con- 
fined to AppendixlAl SectionlVlcontains our conclusions. 

II. MEAN FIELD APPROXIMATION 

The mean-field approach to the Bose-Hubbard model 
introduced by Sheshadri and co-workers El relies on the 
standard approximation 

a 3 a h ~ ( a i) a h + a i( a h) - ( a j)i a h) ( 2 ) 

where aj = (aj) = (a^) is the so-called superfluid pa- 
rameter [32j, to be determined self-consistently. Indeed, 
equation @ allows to recast Hamiltonian Q as the sum 
of terms containing on-site operators only: 

M 

n = (3) 

Hj = —nj{rij - 1) - (/i - Vj)rij 

-t (aj + a+ - a 3 ) ^ r jh a h (4) 

where M is the number of lattice sites. Note that, unlike 
Hamiltonian the mean-field Hamiltonian Ti features 
single boson terms, and therefore it does not conserve the 
total number of bosons. 

A qualitative zero-temperature phase diagram of the 
BH model can be obtained by evaluating the expectation 
value (•) on the ground state of Hamiltonian JSJ. Such 
evaluation must be performed self-consistently, since the 
ground state of TL itself depends on the set of superfluid 
parameters {aj}. In the particular case of homogeneous 
lattices, translational invariance yields aj — a, and one 
is left with (M identical copies of) a single-site problem. 
The resulting phase diagram consists of a superfluid re- 
gion, where a > 0, and a series of Mott-insulator lobes, 
where a — and the local density (nj) is pinned to an in- 
teger value (and hence the system is in an incompressible 
state, d^(nj) = ). The boundaries of these Mott lobes 
have been determined numerically in the original paper 
|17| . while their analytical expression has been reported 
in a quite recent work |2C|. 
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More in general this mean-field approach has been 
adopted for studying the superfluid-insulator transition 
in some inhomogeneous situations. An harmonic confin- 
ing potential Vj oc (j — j ) 2 is considered in Ref. [3^, 
whereas the effect of topological inhomogeneity is ad- 
dressed in Ref. 

Here we are interested in the thermodynamics of the 
system, and we adopt the standard grand-canonical sta- 
tistical mechanic approach, 



(O) 



Tr(e-P H ) 



(5) 



where the trace is evaluated on the whole Fock space. 
Exploiting the site-decoupling of the mean-field Hamil- 
tonian H, Eq. J3J| , the self-consistency conditions become 
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where the traces in the second line are evaluated on the 
single site Fock space. Note that the superfluid param- 
eters ctj can be safely assumed to be real since both Hj 
and cij are real operators. 

In the following sections we illustrate how it is possible 
to determine analytically the critical condition for super- 
fluidity, i.e. for the existence of a stable solution of the 
self-consistency equations © with oij ^ 0. 
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FIG. 1: Finite temperature mean-field phase diagram for the 
homogeneous one-dimensional lattice. The dotted, solid and 
dashed lines are the critical curves, Eq. (|HJ, for T/U = 0, 
T/U = 0.04 and T/U = 0.1, respectively. In every case, the 
superfluid domain is above the relevant critical curve. For 
T/U = 0, the region below the (dotted) critical curve consists 
of two disjoint Mott lobes. For T/U = 0.04 (solid curve) a 
sharp crossover between a normal fluid phase and a insulator- 
like phase is present. The gray areas are the insulator-like 
regions as evaluated setting e = 10 -3 in Eq. 11211 . As discussed 
in the text, the particle density is very close to an integer value 
(also shown) inside these regions. For T/U = 0.1 (dashed 
curve) no insulator-like regions are present (for the same value 
of e). 



III. HOMOGENEOUS CASE 

When Vj = and Tjk = t the system is translation- 
ally invariant and, similar to the above-recalled zero- 
temperature case, Eqs. © reduce to (M identical copies 
of) a single consistency equation. Dropping the site- 
labeling subscripts one gets aj = a = (a) and 

Hj = H = — n(n - 1) - fin - 2dta(a + <z+) (7) 

where d is the dimension of the lattice and we discarded 
the constant term 2tda 2 since it can be factored out from 
both the numerator and the denominator of Eq. (JjJJ. 

In this simple case the self-consistency constraint, 
Eq. © , depends on the (site-independent) superfluid pa- 
rameter a, and it is met when the latter is a stable fixed 
point. It is easy to check that a = is a fixed point of 
Eq. © whose stability depends on the parameters U, /i 
and t. In particular, when a — is unstable, a stable so- 
lution a > is expected and the system is in a superfluid 
state. 

In Appendix^we show that the critical curve defining 
the superfluid domain border is 
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where 
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Note that when a = the expectation value of the par- 
ticle density does not depend on t, being simply 

m^0) = (n)= ^ (10) 

Hence, unlike the zero-temperature case, the a — re- 
gion in general does not yield integer particle density. 
However it can be shown that, for sufficiently low tem- 
peratures, there are \i intervals where p(U, fi, (3) is prac- 
tically constant. This happens when a single term of the 
sums in (fTTTf) . whose label we denote k* , outweighs the 
remaining terms, i.e. when 



e 2 
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where e is a small parameter. The last equation identifies 
the interval 



/,* i ln ( £ ) , t c k* i ln(e) 
k 1 0U < U < + f3U 



(12) 
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FIG. 2: (Color online) Three dimensional representation of the finite-temperature mean-field phase diagram for a homogeneous 
one-dimensional lattice. Dark gray (green): critical surface, Eq. ®. Light gray (yellow): zero-temperature Mott lobes. White- 
rimmed transparent surfaces: boundaries between the insulator-like (below) and normal fluid (above) regions. Figure is 
obtained by cross-sectioning this figure at the relevant values of T/U. 



where p{U, p,fl) ~ k* . More precisely, expanding 
Eq. I|10(l with respect to e and considering only the lin- 
ear order, one gets \p — k*\ « e. In the normal fluid 
region between two subsequent intervals of the form 112|) 
the particle density can be described considering only the 
terms k* and k* + 1 in Eq. (tlTj|) 



p(U,fi,p) 



k* + (k* + l) e 0(M-tffc*) 
1 + e Pbj.-uk*) 



Last equation clearly shows the crossover between two 
subsequent insulator-like regions. Of course, when (3 < 
— 21n(e)/[7 there is no p satisfying the Eq. I|12l) . and the 
plateau-like behaviour of p(U, p, (3) disappears. 

The same line of reasoning allows to obtain quite 
straightforwardly the exact phase diagram of the model 
in the zero temperature limit [20j. Indeed, when f3 — > oo 
only the terms labeled by k — k* survive in Eq. (JSJ), so 
that the critical curve in the p/U-t/U phase diagram is 
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and, according to Eq. 1|12(1 . p G [k* — 1, k*]. 

The boundary of the superfluid domain at different 
temperatures, along with the crossover between the nor- 
mal fluid and the Mott insulator (when present) is shown 
in Figs. Hand [21 for a homogeneous lattice with d = 1. 
As it is evident from Eq. ©, the results for d > 1 are 
obtained by a suitable rescaling of t c . 



IV. SUPERLATTICES 



We now turn to the case of superlattices, where the pa- 
rameters appearing in Hamiltonian are periodic func- 
tions of the site label. Our approach can be applied to 
a generic d-dimensional superlattice, but, for the sake of 
clarity, here we focus on the one-dimensional ^-periodic 
case, Vj = v j+ £ and r jth = {5 h .j + i+8 h , j ^i)T. J+l j l+ (. Note 
that this choice is not merely dictated by the ensuing no- 
tational simplification, but also experimentally relevant 



Since the superfluid parameters mirror the £- 



periodicity of the superlattice, otj 



the self- 



consistency conditions © reduce to £ independent equa- 
tions. As in the homogeneous case, the choice ah = 
for all h's is a fixed point of Eq. ©, and only when it 
is unstable the system is expected to be in a superfluid 
state. According to a standard approach, the stability of 
such fixed point can be discussed based on the spectrum 
of the matrix linearizing the map defined by Eq. in 
the vicinity of the configuration = 0. More in de- 
tail, the fixed point is stable only if the modulus of the 
maximal eigenvalue of such matrix is lower than one. By 
adopting a calculation technique similar to the one of the 
homogeneous case (see Appendix |A"|) . the linearized map 
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turns out to be 
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where, introducing fi h = fi-v h and T hth > = T hh >(5 h h>+i- 

Shh'-l) + T~i,i+l{$hlSh' £ + 5h> lShl), 
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The function Qk(U,n) appearing in Eq. (|15fl is exactly 
the same as defined in Eq. ©. Since F^y is a real 
and positive matrix, Perron-Frobenius theorem j3^ en- 
sures that its maximal eigenvalue <f>M (/?> U, fj) is real and 
positive. Hence, the fixed point ah = is unstable - 
and the system behaves like a superfluid — only when 
t > t c ({3, U,fi) — 4>^. Note that the critical curve sepa- 
rating the superfluid and the normal domains can also be 
defined as the lowest positive t c such that Pp{t~ l ) = 0, 
where -Pf(A) is the characteristic polynomial of matrix 
Fh,h'- 

In the normal phase [t < t c and oth = 0, V h) the local 
density of particles at site h is: 



Ph(U,n,/3) = (n h ) 



J2 k e 0(^hk-^k(k-l)) 



(16) 



Analogously to the homogeneous case, for sufficiently low 
temperatures there exist intervals of p where the local 
particle density is arbitrarily close to an integer value. 
In detail, introducing a small parameter e, \ph{U, (3) — 
k*\ < e if /j, e Mh(k*, e), where 
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and, for any positive integer k*, 
M h (k*,e) 
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This in particular means that the average filling 
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remains very close to the rational value 
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as long as the chemical potential belongs to the interval 



M({k* h },{v h },e) = f]M h (k* h ,e) 



(20) 
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FIG. 3: Schematic representation of the ^-periodic optical 
superlattices explicitly considered in Sec. lIVI The lattice sites 
correspond to the local minima of the optical potential V(x) . 
Upper panel: 1 = 2, Vj = v J+2 Vj+i, r 3 -,j+i = T j+2 ,j+3 ^ 



Tj + l,j+2- 
Tj + l,j+2 



Lower panel: £ 



v J+2 Vj+1, Tj, j + 1 



j+2,j+3, 



j = 3k + 1, k e N. 



where {kt}f l — 1 is a set of non-negative integers such that 
M. ^ 0. If, conversely, oth = but // M., the average 
filling is not close to a rational number and it significantly 
varies with varying fi. That is to say, the compressibility 
is significantly different from zero, and the system be- 
haves like a normal fluid, owing to thermal fluctuations. 

In the zero-temperature limit the normal phase be- 
haviour disappears and a (Mott)insulator-superfluid 
transition is recovered. Similar to the homogeneous case, 
the n/U-t/U zero-temperature phase diagram consists of 
a superfluid domain and a series of insulating Mott-lobes. 
However, the average filling and compressibility within 
these lobes are exactly k* and zero, respectively. There- 
fore, as it is expected Q, we obtain that superlattices 
can display an insulating-Mott behaviour even for some 
critical rational fillings. 

Let us now analyze explicitly the two simple superlat- 
tices generated by the trapping potential schematically 
represented in Fig.|3 The upper panel of this figure cor- 
responds to the simplest choice, namely a superlattice of 
periodicity 1 = 2. In this case the maximal eigenvalue of 
matrix F can be evaluated analytically, and the resulting 
critical value of t turns out to be 



t c (U,n,0) 



F19F2. 



(21) 



h=l 



The ensuing mean-field phase-diagram for a particular 
choice of the parameters is displayed in Fig. 0] As 
mentioned above, rational (actually half-integer) filling 
Mott lobes appear in the zero-temperature phase dia- 
gram (dotted curves). As the temperature increases, the 
regions where the system is in a quasi insulating state 



6 



TIU=0 
771/= 0.02 

7717=0.04 
- - 7717=0.1 




0.1 - 



T/U=0 
T/U= 0.02 
7717= 0.05 
7717= 0.1 



0.05- 




0.5 



1 1.5 

\llU 



FIG. 4: Phase diagram for the superlattice of periodicity I = 
2 with parameters V2k/U = 0.35, V2k+i/U = 0. Note that the 
possibly different values of parameters 77, are conveniently 
absorbed into t c by a multiplicative rescaling. The filled areas 
correspond to the insulator-like phase for T — 0.02 (light 
gray) and T — 0.04 (dark gray) as evaluated setting e = 10 -3 
in Eq. I1201 . The rational numbers denote the particle density 
within the Mott-like domains. 



shrink and eventually disappear, according to Eq. (|20[l . 
It is interesting to observe that these regions may disap- 
pear at different temperatures, depending on their filling. 
This is clearly shown in Fig. where the insulating re- 
gions relevant to T = 0.02 (light gray) and T = 0.04 
(dark gray) are shown. Note that in the latter case there 
are only integer-filling (quasi) insulating regions. 

Matrix F can be analytically diagonalized with a lim- 
ited effort also for the special case of 3-periodic lattice 
illustrated in the lower panel of Fig. [3J where vi = V3 
and 7-1,2 = r 2,3, so that i*i, 2 = F 3i2 , i*i,3 = i*3,i and 
F 2 .i = F 2 ,3. After a straightforward calculation one gets 

t c (U, n, f3) = 2 (22) 

-Fl,3 + J F 1 3 + 8^1,2^2,1 

The relevant phase diagram is shown in Fig. El Note that, 
similar to the previous example, the zero-temperature 
Mott lobes can be divided into two classes depending on 
the relevant particle density, which can be either k or 
k + 1/3, where k 6 N. This simple behaviour is a conse- 
quence of our choice for the local potentials. Of course, 
more structured choices result into a quite richer phase 
diagram, where the particle density in the insulator-like 
regions is an integer multiple of i~ x . Note that some of 
these multiples are excluded if the energy offset between 
any two sites j and h within the same supercell is an inte- 
ger multiple of U. Indeed, in this situation, the intervals 
relevant to sites j and h defined by Eq. (|17|l overlap ex- 
actly. Hence for some set of integers the intersection de- 



FIG. 5: Phase diagram for the superlattice of periodicity 
I = 3 with parameters v 3k /U = v ik+1 /U = 0.35, vsk+2/U = 
0, T3fc,3fc+i = 1/2 and TSfc+i,3fc+2 = T3fc+2,3fc+3 = 1. The 
rational numbers denote the particle density within the zero- 
temperature Mott insulator domains. 



fined by Eq. (|20|l is empty. This is exactly what happens 
in Fig. El where v\ — V3 = 0. Interestingly, it is possible 
to devise superlattices where only fractional critical fill- 
ings are present, provided that the energy offset between 
two lattice sites is larger than U. This is shown in Fig El 
displaying the zero-temperature phase diagram for a su- 
perlattice of periodicity £ = 4 as obtained by numerical 
diagonalization of matrix F, Eq. ©. 



V. CONCLUSIONS 

In this paper we extend the mean-field approximation 
to the Bose-Hubbard model introduced in Ref. ^} to 
include thermal fluctuations. An analytical study of the 
ensuing self-consistency equations allows us to determine 
the exact form of the boundary of the superfluid region 
at any finite temperature. We also quantify the crossover 
from the normal fluid to the insulator-like phase. 

Other than to the homogeneous d-dimcnsional lattice 
considered in the original reference, we apply our method 
to a generic one-dimensional /^-periodic superlattice, giv- 
ing explicit results for I = 2 and a special case of I — 3. 
Results for more complex one-dimensional superlattices 
involve the evaluation of the maximal eigenvalue of a I x i 
matrix. Of course, this must be accomplished numer- 
ically even for relatively small matrix sizes. Neverthe- 
less this approach much less demanding and more pre- 
cise than the (equivalent) fully numerical solution of the 
self-consistency equation. The latter actually involves an 
iterative self-consistent diagonalization for each point of 
the the mesh grid describing the phase diagram. Our 
technique can be extended also to generic d-dimensional 
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FIG. 6: Zero-temperature phase diagram for the superlattice 
of periodicity £ = 4 with parameters v^jU — 1.9, V4k+2/U = 
0.3, Vik+i/U = 1.3, v^k+a/U = 0.0,r fcifc+ i = 1. The rational 
numbers denote the particle density within the Mott insulator 
domains. Note that, as we discussed in the text, for this 
particular choice of the well depths there are no integer critical 
fillings. 



superlattices, where the size of the matrix to be diag- 
onalizcd is s x s, s being the number of sites within a 
supercell. We remark that the boundaries of the super- 
fluid region as evaluated with our method are valid also 
in the zero-temperature limit, thus providing the phase 
diagram for the superfluid-insulator quantum transition. 
In particular this allows to find the analytic expressions 
of the (mean-field) Mott lobes in the case explicitly con- 
sidered above. 
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APPENDIX A 

In this section the analytic expression for the boundary 
of the superfluid domain is derived in detail for the simple 
case of the homogeneous lattice. The generalization to 
the superlattice case is briefly discussed. 

As we mention above, in the homogeneous d- 
dimensional case, the self-consistency constraints, 
Eqs. ijfjj, reduce to a single equation in the variable a, 
i.e. the site-independent superfluid parameter. It proves 



useful to recast such equation in terms of the quantity 

7 = at as 



7 = i- 
where 



Tr ( a e-F* 



Z 



4^ l0g( ^ /(7) (M) 



Z = Tr 



-m 



(A2) 



is the grand-canonical partition function of the single site 
problem. The additional 1/2 factor in Eq. IjAlfl ensues 
from the equality a = (a) = (a + ) G K. 

Let us now prove Eq. (JSJ) by discussing the stability 
character of the fixed point 7 = of the map IjAljl . To 
this aim we truncate the on-site Fock basis considering 
states up to a given number of particles n and denote 
7i n the relevant Hamiltonian matrix. The final result is 
obtained letting n to infinity. 

Introducing the set of eigenvalues of H n , {Ek(l)}k =0 > 
Eq. HAlf) becomes 
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1 ^fc=0 6 ±f_ /AON 

Now, since p n (Ek(^f);^) = 0, where p n (A;7) is the char- 
acteristic polynomial of W„ 



= 



<fpn(-Efc(7)»7) 

dj 
dE k (~?) 



d\p n (A, 7) + dyp n (A, 7) 



\=E k (-f) 



(A4) 



so that it is possible to write 

dE k (j) _ <9 7 p n (A, 7 ) 



(hi 



OxPniKl) 



(A5) 



X=E k ( 1 ) 



Denoting Pn k ^(X;j) the characteristic polynomial of the 
matrix obtained by discarding from Ti. n the rows and 
columns labeled by the set of indices {fc}, and making 
use of the formula for the derivative of a determinant 
|3o| . one gets 



fc=0 



pl fc) (A;7) 



<9 7 Pn(A, 7 ) = 78 d 2 



^fc J3 ( fe ' fc - 1 )(A; 7 ) + P(7) 



Lk=l 



where the polynomial P{"i) is homogeneous, so that 
P(0) = 0. Hence Eq. l|A3jl becomes 



7 = ^7 



where 



5n(A;7) 



1 d T p ra (A,7) 
Ad d\p n (\,j) 



/(7) 



(A6) 



2d 



(A; 7 )+P(7) 



(A7) 
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According to standard treatment, 7 = is a stable solu- one gets 
tion of Eq. only if |<9 7 /(7)| 7=0 < 1, i-e. 



t < 



J2k=o e 



-PE k {0) 



Observing that 



^' } (A;0) = J] (E h (0)-X) 



(A8) 



(A9) 



q n (E k (0):0) = 2d hp^HEkmO) + (* + l)p<£ +1 ' k) (Ek(0); 0)1 (25* (0); 0) 



(A10) 



where we set p 



(n,n+l) 



(A; 7) = Pn (A;7) = 0. Now, 
recalling that Ek(0) — ^k{k — 1) — fj,k one gets 
limn-Kx) q n (Ek(0);0) = Qk(U,p), where the function Q k 
is defined in Eq. (|5J. Therefore, the limit n — > 00 of 
Eq. HA8|> gives the desired result, Eq. 

In the case of superlattices the parameters appearing in 
Hamiltonian are periodic functions of the site labels. 
Here we focus on the one-dimensional ^-periodic case, 
Vj = v 3+i and r jih = (§h,j+i + 8h,j-i)Tj+t,h+l- Since the 
superfluid parameters mirror the ^-periodicity of the su- 
perlattice, otj = aj+e, the self-consistency conditions © 
reduce to t independent equations. Introducing the pa- 
rameters 7/j = ta h (h = l, 7^+1 = 71 and 70 = 7$), 
equations © can be recast as: 



7/i = 1-3 



t 



1 rflog(Z^) 1 d\og(Z h ) 



MM) 



Vi.h+i dj h+1 



where 



Zh = Tr 



(A12) 



and 



H h = -^n(n - 1) - (/i - v h )n - 

(rh,h-ilh-i + T h M+ilh+i)(a + a + ) (A13) 



A procedure similar to that detailedly illustrated in 
the case of homogeneous lattices allows to linearize 



(All) Eq. ETT|) . obtaining Eq. JT2J. 
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